Slave boson theory for transport through magnetic molecules with vibronic states 
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We study the electron transport through a magnetic molecular transistor in the Kendo limit using the slave 
boson technique. We include the electron-phonon coupling and analyze the cases where the spin of the molecule 
is either 5 = 1/2 or 5 = 1. We use the Schrieffer- Wolff transformation to write down a low energy Hamiltonian 
for the system. In the presence of electron-phonon coupling, and for 5=1, the resulting Kondo Hamiltonian 
has two active channels. At low temperature, these two charmels interfere destructively, leading to a zero 
conductance. 

PACS numbers: 72.15.Qm, 73.22.-f 



I. INTRODUCTION 

The technological advances of the last decades triggered 
a systematic study of electronic transport in nanoscale sys- 
tems weakly connected to external electrodes. Confinement 
of electrons at the nanoscale leads naturally to energy level 
quantization and charging effects. In small quantum dots 
Coulomb blockade becomes a dominant effect as shown in 
conductance experiments.'-^ This effect, together with the en- 
hancement of the low-temperature conductance in the valley 
between Coulomb blockade peaks — a signature of the well 
known Kondo effect^'* — makes evident that electron-electron 
coiTelations play a central role in these systems. Since the 
early experiments in quantum dots^M the Kondo effect has 
also been observed in a variety of molecules. The conductance 
of molecular systems like Ceo molecules, Co and Cu com- 
plexes connected to metallic leads or carbon nanotubes, shows 
the signature of Kondo-like behavior . ^'^'^'''^ In these molecu- 
lar junctions, the molecule is coupled to the source and drain 
leads providing a path for charge transport. In some cases a 
third gate electrode is used to create an external potential that 
modifies the energy of the molecular orbitals. The use of nor- 
mal, ferromagnetic and superconducting electrodes gives rise 
to a rich behavior and opens new alternatives for the study of 
the interplay between Kondo screening and electronic coiTela- 
tions in the leads. In addition, transport through highly struc- 
tured molecules like the Mni2, a single molecular magnet, has 
been reported' ' and the Kondo effect has been predicted to oc- 
cur in this type of systems for both half-integer— and integer 
values of the molecular spin^i^ 

Despite of the substantial theoretical and experimental ac- 
tivity devoted to the study of molecular transistors there are 
still many open questions. While the theoretical frame for 
the study of the Kondo effect in semiconducting quantum 
dots and molecular systems is basically the same, there are 
some fundamental differences between them. Namely, vi- 
brational modes can play a central role in the single elec- 
tron transfer through molecules. In fact, molecules dis- 
tort upon the addition or the removal of electrons leading 
to a large electron-phonon interaction. As the Coulomb 
charging energies in these systems can be considerably 
reduced by screening due to the electrodes^ii electronic 



and vibrational energies can become of the same order of 
magnitude generating scenarios where novel effects may 

emerge ,15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32 

The effect of the coupling between electronic excitations 
and vibronic states in molecular transistors depends on the 
symmetry and frequency of the vibrating mode and on the 
strength of the coupling. Symmetric modes with a Holstein 
coupling between quantized vibrations and electronic levels 
may strongly renormalize the molecular parameters reducing 
charging energies'^ and producing anomalous behavior of the 
Kondo temperature versus applied gate voltages. Asymmet- 
ric modes coupled to the tunnel-barriers parameters, like the 
center of mass motion mode, can dynamically open new chan- 
nels for electron transport.^^ Due to the large variety of mag- 
netic molecules that could be incorporated in molecular cir- 
cuits with different chemical environments, it is important to 
characterize and understand the behavior of molecular transis- 
tors with different vibronic and electronic structures. 

In this work we discuss the electronic transport through a 
magnetic molecule with spin 5=1/2 and 5 = 1 in the Kondo 
limit. We use slave boson techniques and include the electron- 
phonon interaction with different coupling constants. We 
show that in the presence of vibrational modes, the low tem- 
perature behavior in the case of 5 = 1 corresponds to a two 
channel Kondo problem. In particular, the zero temperature 
conductance goes to zero. 

The paper is organized as follows: in Sec. II we introduce 
a model Hamiltonian for the system and use the Schieffer- 
Wolff transformation to obtain the low energy Hamiltonian of 
the molecular system for both 5= 1/2 and 5=1. In Sec. Ill we 
use slave boson techniques to find the ground state properties 
of the system in both cases and calculate the conductance. We 
summarize in Sec. IV. 



II. THE MODEL 

We first present the case of a spin 5=1/2 molecule de- 
scribed by the Holstein- Anderson model. The Hamiltonian 
reads i/ = i/o + ^4-ph + ^^hyb, where Hq describes the molec- 
ular degrees of freedom and the electronic excitations of the 
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FIG. 1: Scheme of a molecular device. 
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where creates an electron in a molecular orbital with en- 

,t 

aka 



ergy Em, U is the intramolecular Coulomb repulsion, c 
creates an electron in the k mode of the a lead where a=L,R 
stands for the left and right lead, respectively. The last term 
in Eq. ([TJ describes the molecular vibronic mode of frequency 
COt). The term //e-ph in the Hamiltonian describes the Holstein 
coupling between the electrons in the molecular orbital and 
the molecular vibrations. 
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Finally, the leads-molecule hybridization term is 

-^hyb = -Ttt Va(dlcaka + cl^a^a) 
a.k.a ^ ' 



where is the number of band states. For some special re- 
lations between the parameters of the Hamiltonian, the model 
reduces to known and, in some cases, well studied models: 

i j Vr = Vi, A 7^ and gR=gL = corresponds to a molecule 
with a Holstein mode in a symmetric environment. '^-'^'^^ 

ii) Vr ~Vl, a 7^ with gj^ — gi ^ is a generalized 
Anderson-Holstein model in which vibrations modify the 
molecular energies and the tunneling barriers. 

Hi) Vr = Vl, gR = —gL and A = corresponds to a molecule 
with inversion symmetry and a center of mass motion. ^^'^^ 

iv) Vr = Vl, and gR = —gi with A 7^ describes a molecule 
with no inversion symmetry and a center of mass mode.^^^^ 

In the Kondo regime {—U < Em < 0), when the number of 
electrons in the molecule is well defined, the charge excita- 
tions can be integrated out by means of a Schrieffer-Wolff 
transformation!^ The resulting Hamiltonian is 



H = Y.^kWlkWak+HK 



(4) 



with 



a.p=R,L k.k' 



(5) 



Here vL = i^l kv'^aki') ^ spinor corresponding to the a 
lead, S is the spin operator associated with the molecular or- 
bital and (T are the Pauli matrices. The coupling constants are 
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where Yon = {l/Vn[){X/hcoo)"e-(^/''''^^'/\ Yin = 7o«(« - 
(A/7za)o)^)^tOo/A, EM^£M + ^^/fi(Oo, U = U -IX^/hcoo and 
JiR = Jrl. The electron phonon interaction renormalizes the 
molecular parameters Em and U : the renormalized energy Em 
increases with A while the intramolecular interaction U de- 
creases. This naturally leads to an increase of the Kondo tem- 
perature while it reduces the range of gate voltages where the 
charge state with one electron is stableJ^ For large enough 
A, U becomes negative, which corresponds to an electron- 
phonon induced attractive interaction. In what follows we 
consider the case U >Q.ln this work, we neglect the potential 
scattering term generated by the Schrieffer-Wolff transforma- 
tion. 

It is convenient to make a rotation in the [L,R] space to 
eliminate cross terms. We then define the even (e) and odd (o) 
operators, 

Ceka — UCRi^a + VCLka 

Coka = VCRka-UCLka (8) 

with u^ = {\ + x/l~F)/2, v2 = (1 - VT~^)/2 and F = 



^JlJ [{JRR-JLLf+^JlL 

Hamiltonian reads 
Hk 



In terms of these operators the 

(9) 



ri=e.o k,k' 



where 



Jr^^l;[jRR + Jll ± V {Jrr - JllY + 47; 



'RL 



(10) 



k.a 



with the + and — signs corresponding to the e-channel and 
o-channel, respectively. The Hamiltonian ^ corresponds to 
a two channel Kondo model. In our case the two channels 
character is a consequence of the dynamical nature of the hy- 
bridization operator and the lack of right and left symmetry: 
if gR=gL then J^^^^JrrJll and = . For gR ^ gL the sys- 
tem has always two channels coupled to the molecule's spin. 
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However, as in general Je > Jo, the resulting low temperature 
behavior is dominated by the e-channel. In the next section 
we discuss this behavior in some detail. 

More interesting is the case of a S= 1 molecule. Such a case 
can be described by including an additional spin Sq = 1/2 
that couples ferromagnetically to the spin of the electrons in 
the hybridized molecular orbital described by the first terms 
in Eq. (1). Then, the Hamiltonian includes a term ~Jh Sq ■ S. 
The total spin of the molecule is either zero or one and for 
large values of the Hund-like coupling the singlet state can 
be neglected. Projecting the system onto the total spin S = 
1 state, the Schrieffer- Wolff transformation gives again the 
Hamiltonian of Eq.®, where now the molecular spin operator 
corresponds to 5 = 1 and the coupling constants are half the 
expressions obtained in the previous case with £m replaced 
with eM-3y///4. 

In the context of semiconducting quantum dots, it has been 
shown that the zero temperature conductance for a 5 = 1 dot 
coupled with two channels is G = Gosin^(5i— ^2) were 5\ and 
^2 are the scattering phase shifts at the Fermi level for the two 
channels. Then, we expect molecules with 5=1/2 and 5=1 
to behave quite differently if phonons open the o-channel. Be- 
low we present a slave boson mean-field approach to describe 
this situation. 

Before finishing this section we would hke to emphasize 
that in our model there is only one hybridized orbital so that 
the two channels character arises only from the coupling to 
the phonons. This is different from the case described in Ref. 
ifisll where the authors considered two different orbitals cou- 
pled to the leads. Our model then, corresponds to the simplest 
scenario that allows us to discuss the effect of the electron- 
phonon coupling. In general, when more than one molecu- 
lar orbital is coupled to the leads, there would be more than 
two channels involved, which presumably will lead to a more 
complex behavior 

III. SLAVE BOSON MEAN FIELD THEORY 

For the sake of completeness, we start with a brief discus- 
sion of the spin 5=1/2 case. We review the slave boson 
mean field theory and calculate the conductance. Although 
the results for this case are well known, they are presented as 
a guide to the more interesting case of a spin 5=1 molecule. 

A. The 5 = 1 /2 case 

Following standard procedures, we use a fermionic repre- 
sentation for the molecular spin, 

S, = \{.flf^-flfi), S+^,f]f^, S-=flf^. (11) 
In terms of these operator, the Kondo Hamiltonian reads 

^K=^Y. ■^'?Ll(/T/T-/i/l)(c],^tC^^'T-c],^j,c^n) 

" 77=e,o k.k' 

+ ^(/T/icj,yC^i'T +///Tcj,^c^^,^) . (12) 



The mean field approach corresponds to define bosonic fields 
^■qa = 1/ VNT.k'^l^kaf'^ evaluate the physical quantities at 
the mean field level in these fields. The resulting mean field 
Hamiltonian can be put as— 

"k = - L 4^E(A,-a)4c,,,+//.c. (13) 

where (..) indicates the thermal average. The solution must 
satisfy the condition {^/| /t + ///j ) ^ ^ ^^^^ imposed by 
introducing a Lagrange multiplier. The procedure is well jus- 
tified only when the original model is extended for the spin 
variable a to take .yV different values with o/K — > 00. How- 
ever, it is well known that the approximation captures essen- 
tially the correct low temperature physics even for ^ = 2. 
In the absence of an external magnetic field we may take 
(Ajjl) = {^r}i) = ^7)- The self-consistent solutions for the 
bosonic fields are 

/■^ (0 
A^ = -7r,A,7p j^^dcofico)—^-^, (14) 

where /(ft)) is the Fermi function, p = 1 /2D is the leads den- 
sity of states per spin and 

F = + r„ = np [{JAef + [JoKf] ■ (15) 

For Je 7^ /„ these equations do not have a solution with non- 
zero values for both Ag and A^,. At zero temperature we have 




where r\{—r\) corresponds to the largest (smaller) Kondo cou- 
pling Jrj ■ This solution describes the situation in which one 
channel completely decouples from the molecular spin. In the 
renormalization group language the ratio between the largest 
and smaller couplings flows to infinity and the Kondo screen- 
ing is due only to the channel with the largest coupling. Note 
that for the particular case Je = Jo, the overscreening of the 
5=1/2 spin requires a more elaborate approximation.-'^ 

In what follows we assume that Je > Jo- The tempera- 
ture dependence of the Kondo correlation energy, defined as 
Ek = {Hk) = —4-JeAl, is shown in figure 2. At a temperature 
TKe the correlation energy goes to zero in a singular way. This 
is known to be an artifact of the approximation, TKe should be 
interpreted as a crossover temperature between a high temper- 
ature regime where the spin is essential unscreened and a low 
temperature Fermi liquid regime. According to Eq. (fT4l i the 
Kondo temperature is given by the solution of the following 
equation 

J fD I CO 

l = ^p dco- tanh{—). (17) 

2 J-D CO ZlKe 

The linear conductance through the molecule is calculated us- 
ing a Landauer approach. At the mean field level described 
by Hamiltonian ( fT3] l, the problem reduces to a single resonant 
level centered at the Fermi energy. This structure represents 
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FIG. 2: Temperature dependence of the correlation energy (top) and 
the conductance (bottom) in the 5=1/2 case. Parameters: D = 2, 
pj = 0.22 and Vl = Vr 

the Kondo resonance that leads to resonant tunneUng between 
the source and the drain contacts. 

The transmission function is given by T{<X)) — TrGTiG" 
with = u^Fe and Fl = v^Fg, and G" are the retarded and 
the advanced propagators of the /-fermions, respectively. In 
terms of T{co) the conductance is then, 

The last expression was obtained by using the Green functions 
corresponding to a resonant level of width Eg. For a system 
with inversion symmetry, we have Vl — Vr and = 1 . In 
this situation the zero temperature conductance is the quan- 
tum unit of conductance 2e^/h. As the temperature increases 
Fe decreases and goes to zero at the Kondo temperature Tkb- 
The temperature dependence of the conductance is shown in 
figure 2. As akeady mentioned, the singular behavior at Txe 
is to be interpreted as a crossover. In fact, above the Kondo 
temperature the conductance is small but non-zero and can be 
calculated using perturbation theory.— 

B. The S = 1 case 

From the many possible representations of a spin one op- 
erators, we use a two spin one-half fermion representation 
that has proved to give good results when treated at the mean 
field level. Following Refs. Il39ll40ll we define two states 
with quantum numbers i—1,2 and a =t , i . The operators 
f-^ create a fermion in these states. The physical Hilbert 



space corresponds to the subspace with n,- = «,| + «,| = 1 , with 

'^ia=fi'aficr- The spins of the states 1 and 2 are coupled fer- 
romagnetically with a exchange constant Jfj. Within this sub- 
space, the three components of the spin one operator are given 
by 

Sz = nx]n2] - niin2i 
5+ = 11182 +S'^n2 

S- = 5^n2+ni52 , (19) 

here Sf" are the (raising/lowering) spin operators of state /, 
which are defined as in Eq. (fTTT l. 

As mentioned above, the effective Hamiltonian Hk for 
this case is described again by Eq. (|9]), but the definition 
of the coupling constant Jri is obtained by replacing Em 
Em — ^Jh/^ in Eqs. (|6]l and ^ and dividing Eq. ( fTOl i by 2. 
Defining bosonic fields as 




we can write the transverse part of the spin product as 

" k,k' io 

with /= (; + l,mod2). As in the 5=1/2 case, we treat this 
term in the mean field approximation. This gives 

ia ia 

+ {^n6){^'^a)ni-2{X'^^){X'^,){ni) 

(22) 

where a = a. In the absence of external magnetic fields 
these quantities (A^^j) are spin independent. 
The Kondo Hamiltonian is then given by 

I E (^no)^\kafio+H.c. (23) 

where we have taken (n,) — 1 . The total Hamiltonian also 
includes a term of the form ^//^((n, — 1) +C where jij is 
a Lagrange multiplier introduced to preserve the number of 

/-fermions and C is a constant. The quantities ^A^'fj^ are 

calculated self-consistently. The symmetry of the Hamilto- 
nian in the quantum number ; may suggest that they are i- 
independent. However, the most general solution is 

(A^i)=A„, (A2+)=A,."''- (24) 

with At) a real number. For the one channel case, the fase 
Og can be eliminated by a gauge transformation. However, as 
shown below, in the two channel case the fase difference 89 = 
9e — 9o plays an important role. With the notation defined 
above, the Hamiltonian ( |23] ) reads 
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FIG. 3: Temperature dependence of the correlation energy (top) and 
the conductance (bottom) in the S= 1 case. Parameters: D = 2, pjg = 
0.22, p/o = 0.13 
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' 50 ' 

'sin( — ]Kkafoc+H.c. 



(25) 



where /(.< 
and c'. 



(/la + 

,-i6e/2^ 



^/2a)/V2, (/la " e'^^/za)/ V2, 

^oA:a~'= '-oita- ^'^^ 50=0, the two channels are 
mixed with the same fermion. Then, as in the 5=1/2 case, 
one channel decouples and does not contribute to the energy. 
For 56 = n the mean field problem reduces to the case of two 
channels, each one coupled with a resonant state and both con- 
tributing to lower the energy. It can be shown that for the gen- 
eral case with Jg ^ Jg, self-consistent solutions exist only for 
50 =0 or ±n. Without any loss of generality we can take 
de = 0. Then, for 50 = ±7r we have 



Hk 



E 

r\=e,o 



k.a 



f^^a+H.C. 



where the self-consistent parameters An are given by 



rD 

l = -y^p / d(0f{(0] 

J-D 



(O 



(26) 



(27) 



with r,) = Kpq^ and qj^ =y,)A,j. Accordingly, the zero tem- 
perature correlation energy is Ek = —A-JfA^ — AJ„I^. As the 
two coupling constants are different, there are two character- 
istic energy scales given by the two Kondo temperatures Tkc 
and Tko- The coiTelation energy as a function of temperature 
is shown in figure[3] 

The conductance can be easily evaluated within the mean 
field approximation following the same line as in the 5=1/2 




FIG. 4: Transmission function T[(0) as a function of ft) for a 5 = 1 
molecule. For Tko <T < Tjce the system shows the usual one chan- 
nel Kondo resonance (solid line) while for T < T^,, the odd channel 
becomes active and a zero- transmission dip appears at ft) = (dashed 
line). 



case. The transmission function in this case is given by 
T{(o)=Ti-[TlG"TrG''], where the T's and the propagators 
are 2x2 matrices since we now have two active channels cou- 
pled to the 5 = 1 spin. In the even and odd base we have 



Tr = 7ZP 



qeqe 
qoqe 



qeqo 
qoqo 



qeqe 
-qoqe 



-qeqo 
qoqo 



(28) 



Finally, the conductance can be put as 

dfico) 



2^2 
h 



dco 



d(0 



4«2^2 



(0 + iTe (0 + iTo 



(29) 

here F,) = ^pq^ ■ The conductance as a function of tempera- 
ture is shown in figure [3] Initially the conductance increases 
as T decreases as in the usual one-channel case but it starts de- 
creasing towards a zero value as T becomes smaller than Tko- 
This interference effect between the two channels is clearly 
seen as a dip at o) = in the transmission function (Fig. IHi. 

Interference effects in the context of 'single impurity' 
Kondo physics have been extensively discussed in the 
literature ^ 25 ,4 1 ,42,43 ,44 ^^j. ^^^g jj^g origin of the effect is 
clear as discussed in detail by Pustilnik and Glazman.^^ In 
the low temperature Fermi liquid regime, the mean field ap- 
proach presented here allows for a simple interpretation of 
the interference phenomena: Hamiltonian (|25] | with 50 = ;r 
and non-zero A^- and Ag is equivalent to a two lead junction 
connected by two different paths. In this effective model, the 
phase difference of the two paths is n leading to destructive 
interference and zero conductance at T = 0. 



IV. SUMMARY 

We have analyzed the Kondo Hamiltonians for magnetic 
molecules with vibronic states in molecules with spin 5=1/2 
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and 5= 1 within the slave boson mean field theory. We showed 
that in systems with no L-R symmetry in the electron-phonon 
coupling, there are two channels coupled to the molecular 
spin. For the case of 5=1/2 molecules, the weakest cou- 
pled channel plays no role; at the mean field level it is sim- 
ply decoupled. For the S= \ case, the two channels screen 
the molecular spin at different energy scales. To fix ideas we 
summarize the behavior of two simple cases with identical Vl 
and Vr: 

a) A Holstein-like mode with A 7^ and gR = gi- The 
Kondo Hamiltonian includes a single channel and for the 
5=1/2 case the slave boson mean field approximation re- 
produces the known results. The electron-phonon coupling 
renormalizes the coupling constants. At Tk the conductance 
increases to reach the unitary limit at zero temperature. In this 
case it has been shown that for large electron-phonon coupling 
X the Kondo temperature is weakly dependent on gate volt- 
ages due to the anomalous dependence of coupling constants 
Jaa with the molecular orbital energy £m-^'' For the S=l case, 
a single channel can screen only half of the molecular spin. 
The zero temperature conductance also reaches the quantum 
of conductance value 2e^ /h although strictly speaking the sys- 
tem is not a Fermi liquid since there is an unscreened spin 



one-half at the molecular junction. 

b) A molecule with inversion symmetry and a center of 
mass mode, A = and gR = —gi- The resulting Kondo Hamil- 
tonian includes two channels with different coupling con- 
stants. For the 5=1/2 case the usual behavior is reproduced as 
one channel decouples. In this simple model electron-phonon 
coupling is just due to the modulation of the tunneling bar- 
riers or hybridizations Va, in the lowest order the molecular 
energies and Coulomb repulsion are not renormaUzed. For 
the 5=1 case the spin is screened by the two active chan- 
nels with two characteristic energy scales Tkc and Tk,, ■ For 
Tko -^T T^e the conductance approaches the unitary limit 
and decreases for T < Tk,, being zero at zero temperature. In 
the most general case with A 7^ , and gR ^ gi the conduc- 
tance are reduced by the prefactor Au^v^ . 
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